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Abstract: We introduce a new class of jet algorithms designed to return conical jets with 
a variable AR radius. A specific example, in which AR scales as 1/pt, proves particularly 
useful in capturing the kinematic features of a wide variety of hard scattering processes. 
We implement this AR scaling in a sequential recombination algorithm and test it by 
reconstructing resonance masses and kinematic endpoints. These test cases show 10 - 20% 
improvements in signal efficiency compared to fixed AR algorithms. We also comment on 
cuts useful in reducing continuum jet backgrounds.* 



*The jet algorithms we describe have been implemented as plugins to the Fast Jet package [1, 2]. These 
plugins are available from the authors upon request. 
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1. Introduction 

The high energy frontier of particle physics is probed by hadronic collisions where hard 
scattering events often result in colored partons. Because these partons undergo showering 
and hadronization, one cannot go directly from detector measurements to the four-momenta 
of the scattering. Instead, one must cluster hadrons into jets as an approximation to the 
short-distance kinematics. 

As jets provide an essential window onto hard scattering processes, much recent work 
has focused on the procedures used to construct them: jet algorithms. Many fast, in- 
frared/collinear safe jet algorithms have been developed [3, 4, 5, 6, 7, 8, 9]. These algo- 
rithms cluster hadrons through different mechanisms (via cones or sequential recombina- 
tion) and in different orders (hard to soft, soft to hard, or by angle). To date, however, 
all jet algorithms used for hadronic collisions return jets of uniform characteristic size 
AR = y/(Ar?)2 + (A(/))2 on the (r/, (p) plane. 

In this paper, we modify existing jet algorithms to cluster with variable AR cone sizes. 
To see how this could be useful, consider a resonance at rest in the lab frame decaying into 
two partons. To first approximation, the shower of hadrons resulting from each parton will 
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fall in a circular cone of fixed angular size regardless of the orientation of the decay with 
respect to the beam axis. However, a jet algorithm which uses a fixed cone size in (r/, (p) 
will not reflect this behavior, as fixed AR corresponds to variable angular size. As we will 
see, this can be remedied by letting the cone size of a jet vary as 



where pT is the transverse momentum of the jet with respect to the beam axis. 

Remarkably, this simple modification to the jet cone size finds applications beyond 
simple resonance decays, and can be useful in studying more complicated kinematic struc- 
tures. The algorithms we present are trivially infrared/coUinear safe and are appropriate 
for use at hadron colliders since Eq. (1.1) is invariant under boosts along the beam axis. 
Moreover, as we will discuss in Appendix A, the fact that jet radii become larger at lower 
Pt makes these algorithms especially robust against splittings from detector effects and 
different showering approximations. 

This paper is organized as follows. In Section 2, we discuss sequential recombination 
jet algorithms, extensions to variable AR, and the particularly useful case specified in 
Eq. (1.1) which we call "VR". In Section 3, we quantify the improvement in signal effi- 
ciency by using VR algorithms in three different event topologies: single resonance decay, 
multiple resonance decay, and three-body gluino decay. Section 4 contains a discussion of 
background shaping and jet quality cuts. Section 5 contains our conclusions. Discussions 
of effective jet radii and VR algorithm parameters can be found in the appendices. 

2. New Recursive Jet Algorithms 
2.1 Brief Review 

Modern jet algorithms fall into two general categories: cone-based and sequential recombi- 
nation. We focus on the latter recursive algorithms because their infrared/collinear safety 
is easier to prove, although much of the following discussion could be adapted to cone-based 
algorithms. 

Sequential recombination algorithms begin with a set of four-momenta derived from 
detector calorimeter cells and then recursively combine pairs of momenta into jets. To do 
this, they take a list of initial four-momenta and assign each pair (i, j) a "jet-jet distance 
measure" dij, while each individual four-momenta is assigned a "jet-beam distance mea- 
sure" diB- At each step, the smallest entry in the set of dy and diB is identified. If the 
smallest entry is a jet-jet measure, the two four-momenta are combined into one, their prior 
distance measures removed from the list, and a new entry with the sum four- vector is com- 
puted. If the smallest entry is a jet-beam measure, then the corresponding four-momenta is 
"merged with the beam" and set aside. The algorithm proceeds in this manner until either 
all four-momenta are merged with the beam (an inclusive algorithm) or a predetermined 
distance measure dcut is reached (an exclusive algorithm). After clustering, an inclusive 
algorithm returns the four-momenta merged with the beam whose pT is greater than some 
minimum value, while an exclusive algorithm returns the unmerged four-momenta (that 
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is, those whose diB > dcut)- Here we wih focus on inclusive algorithms because they are in 
more widespread use and because not all algorithms have well-defined exclusive modes. ^ 

It is possible to parameterize the most popular sequential recombination algorithms 
for use at hadron colliders via [9] 

dij = min [p^t'I^p^] R%, diB = P^tI RI (2-1) 

where the values of n for particular algorithms are listed in Table 1, Rij is the AR sep- 
aration between the two four-momenta, and Rq is a free parameter that determines the 
characteristic jet size. Roughly, n > clusters soft items first, n = clusters by angle, and 
n < clusters from hard particles outward. We will be particularly interested in algorithms 
with n < 0, because as emphasized in Ref. [9], such recursive jet algorithms act much like 
an idealized cone-based algorithm for sufficiently negative values of n. 



Algorithm 


n 


kr [4, 5] 


1 


Cambridge-Aachen [6, 7] 





Anti-/cT [9] 


-1 



Table 1: Parameterization of popular sequential recombination algorithms according to Eq. (2.1). 

At lepton colliders, Eq. (2.1) is usually modified by replacing pTi with the energy Ei, 
and Rij with the arc length Sij on the {9, cj)) sphere defined by A5 = \/ (A^)^ -|- (sin 9 Ac/))^. 
This AS" measure will be part of the inspiration for the jet algorithms presented in Sec- 
tion 2.3. 

2.2 Variable Ai? Algorithms 

We now generalize Eq. (2.1) so that 

dij = min [p%pl^j\ R%, diB = pl'lResiPTi?, (2.2) 

where Res{PTi) is a dimensionless number interpreted as an effective jet radius for n < 0. 
Since Eq. (2.2) is invariant to boosts along the beam axis, it is appropriate for use at 
hadron colliders. See Appendix A for more detail on the possible choices for iigff and the 
restriction to n < 0. 

To see why i?cff is an effective radius for n < 0, consider the clustering of two four- 
momenta i and j. These will only be clustered together if 

dij < diB, djB- (2.3) 

^For instance, the anti-kr algorithm [9] assigns smaller diB to harder jets, so these would be merged 
with the beam and not identified as jets if the algorithm were run with a dcut- 
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Let PTi > PTj- As discussed in Appendix A, in order for i?efr to robustly define an effective 
radius, we must take diB < djB- Therefore, the four- momenta are clustered if 




RcsiVTi)^ 




< 1, 



(2.4) 



which requires the AR distance between i and j to be within an effective radius -Rcff- 

To our knowledge, all current algorithms at hadron colliders use a constant Rcs- Al- 
though we will only explore one new algorithm here, in which i?eff oc 1/pt, it is possible 
to invent new algorithms tailor-made to distinct processes. We leave such extensions to 
future work, though some guidelines for choosing R^q appear in Appendix A. In principle, 
a momentum-dependent effective jet size could be implemented in a cone-based algorithm, 
though we expect that the process of finding stable cones (as is required in SIS-Cone [8]) 
would be much more complicated. 

2.3 Introducing VR Algorithms 

We now introduce a particular example of the variable AR algorithms described above 
which we will denote as "VR". To motivate this setup, consider a resonance decaying in 
its center-of-mass frame. If the resonance decays to two jets, these jets will be naturally 
described by circles on the (9,(1)) sphere. This is the setup used to describe jets at e~^/e~ 
colliders where the center-of-mass frame is the lab frame. Unfortunately, the partonic 
center-of-mass frame at a hadron collider is not fixed, and to maintain boost invariance 
one must work with circles in {r],(j)). This setup can present a "Goldilocks" problem in 
choosing jet radii. If one is forced to use a fixed AR radius, then the cones in the central 
region will be large in AS, while those in the forward region will be small. Getting the 
radius just right will necessarily involve a tradeoff. 

Remarkably, the simplest application of the general jet algorithms described above can 
remedy this situation, allowing one to use boost invariant jets with a cone size reflecting 
their true size as measured in the resonance rest frame. This new VR algorithm works by 
letting the effective radius of a jet go as 



where p is a dimensionful constant. 

To see why this captures the desired behavior, again consider the decay of a resonance 
into two partons which shower into jets. In the rest frame of the mother particle, the energy 
E of each jet is fixed by the mass of the resonance with E = m^(,s,/2. The resulting jets 
should have roughly the same characteristic angular size AS, therefore, the quantity EAS 
is approximately the same for both jets. As we will now show, fixed EAS is equivalent to 
Eq. (2.5) to first approximation. 

First note that EAS is a quantity that is approximately invariant under both transverse 
and longitudinal boosts. Intuitively, for small angles, the invariant mass of two massless 
four-momenta is equal to the geometric mean of their energies multiplied by their angular 
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separation. Since invariant mass is a boost invariant quantity, so is EAS. More formally, 
for small angles 



where we have used the fact that E ^ px cosh rj (true for small jet mass) and AS" ~ 
AR/ coshr] (see Appendix B). Since px and AR are invariant under boosts along the 
beam axis, so is EAS to first approximation. We emphasize that this is true even for 
boosts transverse to the beam axis; just define a new (rj, </>) coordinate system along the 
boost axis and go through the same procedure to show that EAS is invariant. 

Putting this all together, we want build a jet algorithm that captures the fact that 
the two jets have constant opening angle in their mother's rest frame. Since EAS is 
approximately boost invariant, this is equivalent to forming jets of constant EAS in any 
convenient frame. In particular, in the lab frame we can use Eq. (2.6) to swap EAS for 
PtAR. Therefore, to get jets of constant EAS, we should choose R^s to scale like 1 /px as 
in Eq. (2.5). This defines the VR jet algorithms. 

From this logic, we expect the parameter p to be proportional to the typical jet size 
measure EAS, and thus proportional to the resonance mass mrcs.- For a more detailed 
discussion of the valid parameter range for p, see Appendix B. 

2.4 Event Topologies v^ith VR-symmetry 

The VR scaling of Eq. (2.5) is applicable whenever there is reason to expect all jets in an 
event to have the same EAS in some frame, and we call these events "VR-symmetric" . 
This is certainly the case for a single resonance decay. Less obvious is that this is true for 
longer cascade decays; even if a cascade involves many intermediate states, it will still be 
VR-symmetric as long as all jets come from the decay of resonances with a common mass. 
VR-symmetry can even be satisfied when there is no actual reconstructable resonance. For 
example, the three-body gluino decay g 2j + Xi would satisfy the requirement toward 
the kinematic endpoint. We will discuss these scenarios in more detail in the next section. 

An important example without VR-symmetry is initial state radiation (ISR). Jets 
from ISR do not have a preferred mother rest frame and so the VR jet cone scaling is not 
appropriate. In the case of resonance production plus ISR, the hardest two jets will have 
the VR scaling, but the ISR will not, so in principle, a hybrid VR/fixed-cone algorithm 
could have better jet reconstruction performance. While we will include ISR in our Monte 
Carlo simulations, we will only study the hardest jets in an event for which VR-symmetry 
is expected to apply. 

Now we would like to address a few caveats to the VR derivation in Section 2.3. We 
derived our expression for Rcs in the small cone limit. In practice, one must account 
for corrections in going to finite-sized jets when choosing reasonable jet parameters (see 
Appendix B). Similarly, for low px the algorithm would return pathologically large jets, 
so one is forced to cut off the jet radius at some -Rmax; so that 
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For small enough -Rmax; Eq. (2.7) effectively defines a kind of hybrid VR/fixed-cone algo- 
rithm as needed for events involving ISR, though in this paper we will not try to optimize 
the value of -Rmax- 

As a secondary issue in the VR derivation, we implicitly assumed that jets will remain 
circular in AR as we boost from the mother rest frame to the lab frame. While this is true 
if the boost is only along the beam axis (as is the case for single resonance production), 
for more complicated event topologies, the particular shape of a jet will depend on the 
orientation of the hard process, and will in general be non-circular. Our algorithm forms 
circular jets in {r], (p) so it will not capture the true jet shape in more complicated decays. 
However, unlike a fixed AR algorithm, the VR algorithms do scale the overall jet size 
appropriately under boosts. Since jets are only conical in a statistical sense anyway, we do 
not expect shape distortions to reduce the efficacy of the VR algorithms. 

We can implement Eq. (2.7) in any existing jet algorithm parameterized by Eq. (2.2). 
However, as discussed in Appendix A, it is only meaningful to define Rqs when n < 0. 
Therefore, it does not make sense to apply it to the kx algorithm where n = 1. We will 
therefore combine it with anti-Zcy [9] (denoted AKT) and Cambridge- Aachen [6, 7] (denoted 
CA), and compare the resulting algorithms, AKT-VR and CA-VR, with the original AKT 
and CA algorithms. A summary of the jet definitions used in our study appear in Table 2, 
and a sample lego plot showing the effects of the VR algorithm in reconstructing an event 
can be seen in Figure 1. 



Shorthand 


n 


Rcs 


AKT 


-1 


Ro 


AKT-VR 


-1 


P/PT 


CA 





Ro 


CA-VR 





P/PT 



Table 2: The four jet algorithms used in this study, as parameterized by Eq. (2.2). For the VR 
algorithms, we also impose a maximum jet radius as in Eq. (2.7). 



3. Jet Reconstruction Performance 

Here we use Monte Carlo simulations to compare the AKT and CA algorithms with their 
VR cousins. We investigate three different kinematic scenarios to see how the algorithms 
perform with both large and small jets. We focus only on signal efficiency in this section, 
and turn to background rejection issues in Section 4. 

Our signal and background samples have both been generated in Pythia 6.4. 14 [10], 
with parton- level signal events generated in MadGraph 4.4.5 [11]. We use nominal LHC 
beam parameters (14 TeV proton-proton collisions). Final state hadrons are grouped into 
6r] X 6(p = 0.1 X 0.1 calorimeter cells between —3 < rj < 3 and assigned massless four- 
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Figure 1: The same event reconstructed by anti-fc^ (left) and its VR modification (right). Note 
that in going to the VR algorithm, the high-pT jets (dark blue, green) have been reduced in size 
while softer jets (yellow, purple, light blue) have grown. In this example, only the two harder jets 
are expected to exhibit VR-symmetry, and the softer jets are saturating the i?niax =1-0 constraint. 

momenta based on the calorimeter energy. These calorimeter cells are the starting point 
for the recm'sive jet clustering. 

We use the Fast Jet 2.3.4 [1, 2] package for the AKT and CA algorithms, and we 
wrote new Fast Jet plugins for the AKT-VR and CA-VR algorithms. For each kinematic 
scenario, we scan over a range of jet parameters to optimize the jet algorithm performance. 
To keep the comparison fair, we limit the maximum effective AR of the VR jet cones 
using i?max as in Eq. (2.7), and scan the Rq parameter of the fixed cone algorithms from 
to -Rmax- In the three cases below, we find a universal improvement in using the VR 
algorithms over their fixed AR cousins. 

3.1 Resonance Decays Without Background 

The simplest test of a jet algorithm is resonance reconstruction without standard model 
background. We consider resonances with backgrounds in Section 4. Here we consider 
the scenario of a color-octet scalar X, of negligible width, in the process gg ^ X ^ 
gg? We scan the jet parameters up to a maximum radius iimax = 1-5, and optimize the 
parameters to maximize the percentage of events reconstructed in a narrow mass window 
{nix i 25 GeV) around the true resonance mass.'^ The results of this optimization are 
shown in Table 3 for four different values of mx- 

The resonance invariant mass plots from this analysis can be seen in Figures 2 and 
3. The results indicate a uniform improvement in going from the original algorithms to 
their VR variants: the reconstructed resonances are narrower and taller. The relatively 
large cone sizes found in the optimization should not be troubling as similar results were 
found in Ref. [12]. The optimized choices of p displayed in Table 3 also make intuitive 
sense. We expect most of these jets to have a pT slightly below half the resonance mass, 

^The X couples to gluons via the operator Tr(XG'^^G'"^). 

^The ±25 GeV mass window was chosen by hand to approximate the width of the reconstructed peaks 
after showering and hadronization. It is not related to the perturbative resonance width, which is zero, or 
calorimeter smearing, whose effects we have not included. 
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Figure 2: Invariant mass distributions for the X ^ gg resonance decay using the optimized AKT 
and AKT-VR parameters in Table 3. The distributions have the same normahzation, and the y-axis 
is in arbitrary units (A.U.). The VR algorithms yield a better reconstruction, both in the height 
and width of the resonance. 
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Figure 3: Same as Figure 2 but comparing the CA and CA-VR algorithms. 
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Algorithm 


500 GeV 


1 TeV 


2 TeV 


3 TeV 


AKT ^ AKT-VR 
CA ^ CA-VR 


18% (0.9,200) 
17% (0.9,175) 


14% (1.0,450) 
14% (1.0,400) 


10% (1.2,1000) 
7% (1.2,1000) 


8% (1.3,1500) 
9% (1.2,1500) 



Table 3: Percentage increase in the number of events reconstructed in the mass window mx ± 
25 GeV for the VR variant over the original algorithm. The numbers in parenthesis are the optimized 
parameters for the original and VR variant (i?o f^rid p, with p in GeV) respectively. We see that 
the effective Ai? ~ p/rnx of the VR algorithms is comparable to the fixed Ai?. 



so for p ~ mx/2, Rcs will be 0(1), close to the optimized value found for the fixed Ai? 
algorithms. Note that the values of p are consistent with the discussion in Appendix B. 

It is interesting that the improvements offered by the VR algorithms are greatest for 
small resonance masses. The reason is that signal degradation can come from out-of- 
cone corrections and contamination from the underlying event (including ISR). For larger 
resonance masses with high pT jets, underlying event corrections are proportionally smaller, 
so one need only make the cone size sufficiently large to capture the bulk of the resonance 
signal. For smaller resonance masses with low pT jets, the underlying event corrections are 
more important, and the VR algorithms do a better job balancing the need for large cones 
to capture the signal against the need for suppressing contamination. 

As a side note, the AKT algorithm does perform better than the CA algorithm, but 
by only 1-5% in the cases considered. Thus, the improvement shown by the VR algorithms 
is greater than the difference between the two fixed cone algorithms. 

3.2 Longer Decay Chains 

The scenario considered above involved a simple event topology and relatively large jet 
cones. Moreover, the jets were only boosted along the beam axis, so we expect them to be 
approximately circular in (77, <j)). Here we test a more complex scenario where the final state 
is more crowded (requiring smaller cones), and the final state jets can be boosted along 
a transverse axis. As discussed before, when jets are produced in a transversely boosted 
frame we expect the VR algorithms to capture their scaling, but not their exact shape. 
Here we will see that by accounting only for this scaling we are able to realize significant 
performance improvements. 

We use two color-octet scalars, X and V, again with negligible width, of mass mx = 
3 TeV and my = 500 GeV, decaying to jets via gg ^ X ^ YY — > gggg. After clustering 
the jets, we optimize the algorithm parameters to maximize the number of events where 
two pairs of jets each reconstruct my within 25 GeV of its true value, and all four jets 
reconstruct mx within 50 GeV. Here we limit the maximum cone size to iZmax = 1-0. The 
results are shown in Table 4. We again see a universal improvement in reconstruction. The 
reconstructed distributions for mx and my can be seen in Figures 4 and 5, respectively. 

3.3 Three-Body Gluino Decay 

As we remarked in Section 2.3, we expect the VR algorithms to improve jet reconstruction 
as long as there is some reference frame in which the jets have the same energy and opening 
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Figure 4: Invariant mass distributions of the heavy resonance X for the X YY gggg 
cascade decay scenario using the optimized parameters in Table 4. The distributions have the same 
normahzation, and the ?;-axis is in arbitrary units (A.U.). Tliese plots are made after insisting that 
pairs of jets reconstruct the lighter resonance mass niY within 25 GeV. 




Figure 5: Same as Figure 4 but for the lighter resonance Y. 
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Figure 6: Invariant mass distributions of dijets in three-body gluino decay for the optimized 
parameters in Table 5. The distributions have the same normalization, and the y-axis is in arbitrary 
units (A.U.). For the chosen spectrum we expect to see an endpoint at 900 GeV. It can be seen 
that the VR algorithms fall more sharply near this endpoint than the fixed cone algorithms. 
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Algorithm 


X ^YY 


AKT ^ AKT-VR 
CA ^ CA-VR 


15% (0.7,450) 
23% (0.6,450) 



Table 4: Percentage increase in the number of reconstructed events for the process gg X ^ 
YY gggg- We insist that pairs of jets reconstruct my ± 25 GeV, and that four jets reconstruct 
mx ± 50 GeV. The numbers in parenthesis are the optimized parameters for the original and VR 
variant (Rq and p, with p in GeV) respectively. 

angle. Here, we demonstrate a useful application of this effect with a gluino decaying to 
a neutralino lightest supersymmetric particle (LSP), g ^ q + q + Xi- If intermediate 
squark is on-shell, then we expect no such reference frame, as the "upper" and "lower" 
jets have different preferred mother frames. However, if the squark is off-shell, then the 
process is a three-body decay, and the gluino rest frame is a preferred frame for both jets. 
Moreover, we expect VR-symmetry to be enhanced near the kinematic endpoint where in 
the gluino rest frame, the two jets are back-to-back while the LSP is at rest. Since gluinos 
are typically produced with small transverse boosts in the lab frame, we do not expect 
large distortions of the jet shape. 

To test our algorithm, we consider the associated production of gluinos with a neu- 
tralino LSP via ff^g + Xi and use a spectrum where rrig = 1 TeV and m^o = 100 GeV. 
We choose the associated production channel for demonstration purposes to eliminate the 
combinatorial confusion present in gluino pair production. To ensure a three-body decay, 
we have lifted the squarks out of the spectrum by placing them at 5 TeV. With this spec- 
trum, we expect to see an endpoint in the dijet invariant mass distribution at 900 GeV. 
The more accurate a jet algorithm, the better it can reconstruct this endpoint. Therefore, 
we define the measure of reconstruction performance to be the difference in the number of 
events reconstructed in 100 GeV windows above and below 900 GeV. We have optimized 
this measure with -Rmax = 1.5. The improvement of the VR variants is shown in Table 5 
and Figure 6. 



Algorithm 


9 ^ q + q + Xi 


AKT ^ AKT-VR 
CA ^ CA-VR 


14% (1.2,600) 
7% (1.3,650) 



Table 5: Improvement, measured as the difference in the number of events reconstructed in 
100 GeV bins on either side of the endpoint, in using VR algorithms to reconstruct three-body 
gluino decays. The numbers in parenthesis are the optimized parameters for the original and VR 
variant (i?o and p, with p in GeV) respectively. 
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4. Resonance Decays With Background 

In the previous section, we showed that VR algorithms generically lead to improvements 
in signal reconstruction for events meeting our VR-symmetric criteria. It is natural to 
wonder how the VR algorithms will handle background events. The VR algorithms have 
a dimensionful parameter p, unlike their fixed cone counterparts, so some shaping of the 
background might take place. However, we will see that this is not the case when one 
imposes reasonable jet quality cuts.^ 

To understand the effect of the VR algorithms on background, we consider the color- 
octet resonances previously described in Section 3.1 on top of a background of QCD dijets. 
We optimize our algorithm parameters to maximize signal significance, defined as S/^/B, 
where S and B are the signal and background cross sections, respectively, in a given mass 
window. In addition, we remove the enormous low-p^^ QCD background by requiring the 
two hardest jets each satisfy pT > m^cs./^- 

Let us first consider the dijet background without any kind of quality cut, as shown by 
the solid-lined histograms in Figure 7. Here we see significant background shaping by the 
VR algorithms: they are taking low pT jets from events with small dijet invariant mass, 
assigning the jets large cones, and thus pushing the events to a higher dijet invariant mass 
compared to the fixed cone algorithms. Because the QCD background exhibits a steeply 
falling distribution in invariant mass, this increases the normalization of the background 
in the mass windows shown in Figure 7. However, this problem is not particular to the VR 
algorithms, and fixed cone algorithms with large Rq will inevitably increase the background 
to an analysis by having cones that swallow soft particles toward their edges. 

Such contamination can be effectively controlled by jet quality cuts, which can take 
many different forms. For example, we could place cuts limiting jets to the central region of 
the detector. Another possibility would involve systematically subtracting a four-momenta 
related to the "catchment area" of a jet [13]. Here, we explore another approach by 
requiring the energy- and pT-weighted jet centers approximately coincide. In practice, this 
means imposing 

AR{pE,Ppr) <6, (4.1) 
where Pe and Pp^. are the energy- and py-weighted jet centers defined by 

PE = ^EiPi, PpT = ^VTiPu (4.2) 

i i 

and Pi is the massless four-vector of unit length corresponding to the ith calorimeter cell. 

It is reasonable to impose the jet quality standard of Eq. (4.1) because, as discussed in 
Ref. [14], we expect jets from a massless parton to have a small, hard central core. For our 
analysis we choose 5 = (0.04, 0.025, 0.015, 0.01) for mx = (500 GeV, 1 TeV, 2 TeV, 3 TeV). 
The result of imposing such a cut on the background is shown in Figure 7. These 5 values 
were chosen by estimating the characteristic size of AR(Pe, Ppj,) for each mass window. 
We see that this helps significantly in reducing background. If such a quality cut reduced 

* These quality cuts were not used in our signal-only analysis. 
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Figure 7: Invariant mass distributions of the QCD dijet background for the AKT and AKT-VR 
algorithms, plotted in windows corresponding to the four fiducial resonances. The distributions 
have the same normalization, and the y-axis is in arbitrary units (A.U.). The solid (dotted) lines 
show the background without (with) quality cuts. The red lines show the original algorithm, and 
the black lines show the VR modification. The results are qualitatively similar for CA and CA-VR. 
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Figure 8: Same as Figure 7 but for the resonance signals. 



-13- 



the signal as much as the background it would be of no use, but as can be seen in Figure 8, 
these cuts have a milder effect on the signal than on background, since the signal jets do 
have hard central cores. One expects that the S parameter could be optimized to increase 
the signal yield, but we will not pursue this possibility here. 

Now we employ the quality cuts and present the improvement in S/^fB significance 
by using the VR algorithms. We consider a narrow mass window {mx ± 25 GeV) around 
the resonance as in the signal-only study from Section 3.1. The results are shown in 
Table 6, and we again see a universal improvement in going to VR. The jet parameters are 
comparable to the ones seen in the signal-only study of Table 3. 



Algorithm 


500 GeV 


1 TeV 


2 TeV 


3 TeV 


AKT ^ AKT-VR 
CA ^ CA-VR 


19% (0.9,225) 
10% (0.8,200) 


23% (0.9,450) 
18% (0.9,450) 


18% (1.2,1000) 
15% (1.2,1000) 


21% (1.3,1500) 
18% (1.3,1500) 



Table 6: Percentage increase in S/\/l3 significance for resonance reconstruction over background 
within a, fixed mass window mx±25 GeV. The numbers in parenthesis are the optimized parameters 
for the original and VR variant (i?o and p, with p in GeV) respectively. 

While accurate mass reconstruction is important for measuring the properties of a 
hadronic resonance, for discovery of a resonance one simply wants to increase S/^fB. 
Thus, we also consider a floating mass window, of variable size around mx-, chosen to 
maximize S/ ^/B for each jet parameter. Results for this study are shown in Table 7. Here 
the improvement in going to VR is somewhat less than in the fixed window case, but still 
notable. 



Algorithm 


500 GeV 


1 TeV 


2 TeV 


3 TeV 


AKT AKT-VR 
CA ^ CA-VR 


10% (0.6,225) 
3% (0.6, 200) 


21% (0.9,500) 
16% (0.8,450) 


14% (1.2,1000) 
11% (1.2,1000) 


8% (1.4,1600) 
13% (1.4,1700) 



Table 7: Percentage increase in S/^/B significance for resonance reconstruction over background in 
a variable mass window. The numbers in parenthesis are the optimized parameters for the original 
and VR variant (i?o and p, with p in GeV) respectively. 

5. Conclusion 

We have constructed a new class of jet algorithms in which jet radii Ai? become functions of 
jet pt- Although we have only explored the simplest new algorithm in this class, the rules 
described in Appendix A should provide a well-defined sandbox in which more complicated 
variants can be explored. Surely undiscovered algorithms exist within this framework that 
can further improve jet-based analyses at hadron colliders. 
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In this paper, we focused on the simplest variable AR algorithm, denoted VR, in which 
AR scales as 1/pt- Remarkably, this scaling captures the physical size of jets for many 
diverse processes. To test this algorithm, we developed a sequential-recombination imple- 
mentation using the Fast Jet [1, 2] package, extending the Cambridge- Aachen [6, 7] and 
anti-Zcy [9] algorithms to VR variants. In the analysis of single/multiple resonance decays 
we routinely saw 10 - 20% increases in signal efficiency when using the VR algorithms. 
A similar improvement appears in reconstructing the kinematic endpoint in three-body 
gluino decay. 

In order to use the VR algorithms in the presence of a significant continuum jet back- 
ground, we needed to impose jet quality cuts. We emphasize that these cuts are useful in 
and of themselves, and have strong physical motivation, since jets from hard partons have 
different substructure compared to jets formed from soft radiation. After quality cuts are 
imposed, the VR algorithms outperform their fixed cone cousins by 10 - 20% in statistical 
significance. 

Given our success in developing a jet algorithm for VR-symmetric event topologies, 
one could imagine developing more powerful "designer" jet algorithms to improve signal 
acceptance or background rejection for specific physics scenarios. The original anti-ZcT 
algorithm emphasizes the utility of recursive jet algorithms, since it returns idealized cone 
jets in an infrared/collinear safe manner. The AKT-VR variant emphasizes the flexibility 
of recursive jet algorithms to adapt to different event topologies. One suspects that by 
including additional global event data (or even jet substructure [15, 16, 17, 18, 19] data) 
in the definitions of dij and diB, one could better extract information about specific hard 
scattering processes despite the complicated hadronic environment. 
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A. Effective Jet Radii 

It is important for jet algorithms to be infrared/collinear safe, since this ensures insensitivity 
to detector effects and is necessary for meaningful higher order calculations [20]. Since 
variable AR algorithms are extensions of the jet parameterization from Eq. (2.1), it is easy 
to check that they inherit infrared/collinear safety from standard recursive jet algorithms. 
However, the purpose of variable AR algorithms is to define jets with an effective radius 
RcS, and infrared/collinear safety alone is insufficient to guarantee a reasonable notion of 
jet radius. Therefore, we will impose a stronger condition, requiring our algorithms to be 
"collinear robust". 
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We focus our attention on jet algorithms with n < 0, because these are the only 
algorithms for which the concept of an effective radius is meaningful. The reason for this 
is simple: to talk about an effective radius, the jets must be approximately circular, and 
the only way to achieve this with a sequential recombination algorithm is to start with a 
central core and add on to it. This is achieved for n < because clustering begins with 
the hardest four-momenta.^ The kx algorithm does not satisfy this condition, since the 
algorithm clusters from soft objects to hard ones, and the final shape of a jet is only 
determined in the last few stages of clustering. As a result, kx jets assume non-circular 
shapes, and we are unable to assign these jets a meaningful effective radius. 

It is useful to consider three degrees of coUinear robustness. The weakest of these 
demands that the clustering of two four-momenta is controlled by the effective jet radius of 
the harder one. This ensures that the hardest jets in an event have well-defined jet radii, 
and non-circular effects only appear at lower pT- Consider two four- vectors i and j, with 
PTi > PTj- For these two four- vectors to be clustered together, they must satisfy both 
dij < diB and dij < djs- For n < 0, the first inequality yields 

dij ■'^ij 

diB ResiPTi) 

and Rcs defines an effective radius as desired. To make sure that the second inequality 
does not affect the clustering and ruin the interpretation of an effective radius, we require 
diB < djB, which implies 

P^iResiPTi) <PTjRes{pTj). (A.2) 

This is the requirement for a minimal degree of collinear robustness, and is satisfied by 
both the fixed cone and VR {AR oc 1/pt) algorithms. 

A stronger version of collinear robustness is that the jet algorithm should be insensitive 
the resolution of the calorimeter. Even if Eq. (A.2) is satisfied, it is possible for the jets 
formed with a fine calorimeter resolution to be different from the jets formed with a coarser 
resolution. This can happen, for example, when n = — 1 and ReS oc pT, where the effective 
jet radius can become comparable to the calorimeter resolution, and a coarse calorimeter 
can form a single jet out of two four-momenta that would not otherwise be clustered using 
Rcs- To guard against this pathology, one can impose 

Rcs{Pt) > -Rcalorimoter (A. 3) 

for all values of px- For the VR algorithm, this constraint only applies for very high px 
jets, and in practice is never needed for typical beam energies and calorimeter resolutions. 

The strongest version of collinear robustness requires that a jet algorithm should be 
insensitive to macroscopic splittings within the jet radius. This ensures that reasonable 
rearrangements of a jet's substructure do not cause the jet to be reconstructed differently.^ 

^For the marginal case where n — (so clustering is by angle), this still works because the hard center 
of a jet sees a high concentration of radiation at small angular separation. 

®By "reasonable rearrangement", we have in mind situations where an input four- momenta is broken 
down into smaller pieces within the original effective jet radius. 



(A.l) 



This requires that two four-vectors pi and pj should be clustered together if they lie with 
the effective radius defined by the sum four- vector pi + pj . Again assuming that pxi > PTj , 
the desired condition is that dij < diB whenever Rij < Res{PT{i+j))i which implies a 
"shrinking cone" requirement 



RcsiPTi) > Res{PT{i+j)] 



(A.4) 



This is satisfied by the VR (and fixed cone) algorithms, since Resipr) is monotonically 
decreasing (or constant) as pT increases. Note that Eq. (A.4) implies both Eq. (A. 3) and 
Eq. (A.2). 

We emphasize that these requirements of collinear robustness are not necessary for 
the theoretical consistency of the variable AR algorithms. From the point of view of 
perturbative infrared/collinear safety, there is no singularity associated with finite angle 
splittings. However, without some version of collinear robustness, the effective jet radii 
would be difficult to understand, since the jet algorithm would be overly sensitive to the 
precise four- vectors used in the reconstruction. The shrinking cone requirement of Eq. (A.4) 
is the conceptually simplest way to enforce collinear robustness, since it does not require 
defining an i?caiorimctcr and is independent of the precise choice of n. 



B. Valid VR Parameter Range 

Here we discuss the operational range of the jet algorithms discussed in Section 2.3. We 
are interested in the range of p for which the VR algorithms will correctly reflect the true 
jet shape. Note that the following results are only for guidance in choosing the value of p, 
and not a strict set of rules. 

In deriving the VR algorithms, we made use of the fact that distances in AS* were 
linearly related to distances in AR via 



AR ~ AS cosh r]. 

Recalling the definition of pseudorapidity i] in terms of the polar angle 6 



cosh T] 



1 



sm( 



de 

dr] 



sm ( 



it is easy to see that locally, 

dS'^ = de^ + sin^ 9 dd)"^ 



sin^ e {djf + d(p'^)= sin^ 6 dR^ 



which is equivalent to Eq. (B.l). 

The approximation in Eq. (B.l) is valid as long as 



de 




d^e 


dr] 


> -Ar] 


drf 



which conservatively implies (taking A(j) 0) 



Ai?< 2 



de 




d^e 


dr] 


1 


dr]^ 



cos e 



(B.l) 



(B.2) 



(B.3) 



(B.4) 



(B.5) 
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For the VR variants we are considering with AR = p/pr, distortions in the jet shape in 
going from AS" to AR will be sufficiently small when 



COS y 

As a general rule of thumb then, we expect the VR algorithms to correctly reproduce 
the size of jets as long as p < 2pT for most events reconstructed. Because one must use a 
cutoff on i?eff to remove spuriously large cones, events which violate this condition are not 
necessarily grossly incorrect; they will be reconstructed as if by the CA or AKT algorithms 
working with i2max- 
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